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Abstract. We consider a multi-species generalization of the totally asymmetric 
simple exclusion process (TASEP) with the simple hopping rule: for a and /3th- 
class particles (a < the transition a/3 occurs with a rate independent 

from the values a and /3. P. A. Ferrari and J. Martin (2007) obtained the 
stationary state of this model thanks to a combinatorial algorithm, which was 
subsequently interpreted as a matrix product representation by Evans et al 
(2009). This 'matrix ansatz' shows that the stationary state of the multi-species 
TASEP with N classes of particles (A/"- TASEP) can be constructed algebraically 
by the action of an operator on the {N — 1)-TASEP stationary state. Besides, 
Arita et al. (2009) analyzed the spectral structure of the Markov matrix: they 
showed that the set of eigenvalues of the A/"- TASEP contains those of the (N — 1)- 
TASEP and that the various spectral inclusions can be encoded in a hierarchical 
set-theoretic structure known as the Hasse diagram. Inspired by these works, 
we define nontrivial operators that allow us to construct eigenvectors of the N- 
TASEP by hfting the eigenvectors of the ( - 1)-TASEP. This goal is achieved by 
generalizing the matrix product representation and the Ferrari-Martin algorithm. 
In particular, we show that the matrix ansatz is not only a convenient tool to 
write the stationary state but in fact intertwines Markov matrices of different 
values of N. 
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1. Introduction 

Interacting particle systems are mathematical models used to study collective 
properties of N-body systems evolving with time [45 . Replacing the actual evolution 
rules (given by a classical or a quantum Hamiltonian) by a stochastic dynamics leads 
one to formulate many properties of the system in terms of probability theory and 
allows the use of powerful mathematical methods p4 . Hence, universal properties of 
nonequilibrium statistical physics are successfully explored through the investigation 
of stochastic processes, that are governed by simple dynamical rules at the microscopic 
level but that display a rich macroscopic behavior |41]|43]. When defining such models, 
some specific features of the original physical problem have to be discarded, but if the 
relevant and fundamental characteristics are retained (such as symmetry properties 
and conservation laws) it is expected that the large scale and large time behavior of the 
system will be described correctly [25 . Building a simple representation for complex 
phenomena is common procedure in statistical physics, leading to the emergence of 
paradigmatic models: the harmonic oscillator, the random walker, the Ising magnet 
and so on. These "beautiful models" often display wonderful mathematical structures 

IS HI]. 

In the field of nonequilibrium statistical mechanics, the asymmetric simple 
exclusion process (ASEP) has reached the status of a paradigm [Mj [29l |35]. The 
ASEP is a lattice-gas model of interacting particles. Each of these particles is a 
random walker that hops from a site to one of the neighboring locations, but a move 
is allowed only if the target site is empty. This exclusion constraint mimics short-range 
interactions amongst particles. In order to drive this lattice gas out of equilibrium, 
non-vanishing currents have to be established in the system and this can be achieved 
by various means: by starting with a non-uniform initial condition that takes an 
infinite amount of time to relax, by coupling the system to external reservoirs that 
drive currents through the system (transport of particles, energy, heat, etc.) or by 
introducing some intrinsic bias in the dynamics that favors motion in a privileged 
direction, each particle being an asymmetric random walker that drifts steadily along 
the direction of an external driving force. In particular, the case where moves are 
allowed in only one direction is said totally asymmetric (TASEP). 

The ASEP has been invented several times and in different contexts due to its 
simplicity. It was probably first proposed as a prototype to describe the dynamics of 
ribosomes along RNA [36 . In the mathematical literature, Brownian motions with 
hard-core interactions were first studied by Spitzer [44 who coined the name exclusion 
process. The ASEP also appeared naturally in the description of transport processes in 
systems with strong geometric constraints such as macromolecules transiting through 
capillary vessels and that cannot overtake each other [32 , or anisotropic conductors 
known as solid electrolytes where electrons hop from a vacant location to another 
and repel each other via Coulomb interaction [17. Popular modern applications of 
the exclusion process include molecular motors that transport proteins within the 
cells along actin filaments [30] and, last but not least, the ASEP and its variants are 
ubiquitous in discrete models of traffic ffow [40 . 

From the mathematical point of view, the ASEP is one of the simplest but 
nontrivial models for which the hydrodynamic limit can be rigorously proved. At 
large scales, the distribution of the particles of the ASEP emerges as a density field 
that evolves according to the Burgers equation with a vanishingly small viscosity 
[3TJ |45l [48] . The Burgers equation is the textbook prototype for shock formation: a 



Recursive structures in the multispecies TASEP 



3 



smooth initial distribution can develop a singularity (a discontinuity) in finite time. A 
natural question that arises is whether this shock is an artifact of the hydrodynamic 
limit or if, under some specific conditions, the original ASEP does display some 
singularity at the microscopic scale [H |4l [iTJ [12] [27]. This question was answered 
positively [16] : a shock does exist at the level of the particle system and its width is of 
the order of the lattice size. However, defining precisely the position of the shock at 
the microscopic level requires a trick which is achieved by introducing a new type of 
particle called a second-class particle^ which amounts to coupling two TASEP models 
which initially differ at a single location [33^. This second-class particle, denoted by 
2, has the same dynamics as a normal (or first-class) particle, denoted by 1, but first- 
class particles treat it as a hole (denoted by 0). The local dynamical rules thus take 
the following simple form: 

10^01 12^21 and 20^02, (1) 

where all transitions occur with a same rate. It was proved rigorously in ^20] Hi] that 
the dynamics of the unique second-class particle mimics the motion of the shock in an 
infinite system. Alternatively, by considering a finite density of second-class particles 
in a periodic system, the shape of the microscopic fronts were exactly calculated [16]. 

A straightforward generalization of the two species case of first- and second-class 
particles is to the multispecies exclusion process where there is a hierarchy amongst 
N different species (or classes) of particles [ij |2l [37] (A^-ASEP): a particle of class 
a views particles of classes /3 > a as holes and is viewed as a hole by particles of 
classes /3 < a; in other words, "higher" classes have lower overtaking priority, as 
will be defined precisely in the next section [|j The stationary state of this simplest 
multispecies generalization of the ASEP is highly nontrivial: it was first constructed 
by Pablo Ferrari and James Martin [23, 2i in the TASEP case through a mapping to 
queueing processes. Their construction was inspired by earlier combinatorial results 
[5] |22l [17]. In [iHl [38], the Ferrari-Martin algorithm was re-expressed as a matrix 
product representation in which the weight of a configuration in the A/"- ASEP is written 
as a linear combination of weights of configurations of the {N — 1)-ASEP. This matrix 
(product) ansatz hence allows one to derive steady-state properties of a given model 
knowing those of a simpler system. 

More general spectral properties of the multispecies ASEP on a ring were 
investigated in [6 where models with different total number of species were related to 
each other in a different manner. In that work, the key observation was that particles 
of two consecutive classes cannot be distinguished by particles of other classes. Thus, 
identifying two consecutive classes n and n + 1 defines a natural projection from the N- 
ASEP onto the {N — 1)-ASEP. Because one can identify any two consecutive classes of 
particles, there are N different such projections. These mappings together with their 
commutation relations endow the set of all possible models with a poset structure 
represented by a Hasse diagram [6 . The existence of this structure leads to canonical 
inclusion relations and duality in the spectrum of the Markov matrix. In particular, 
through the identification mapping, the eigenvectors of the A/'-ASEP either vanish or 
project onto eigenvectors of the {N — 1)-ASEP. However, these projections blur some 
essential information (namely, the difference between the two consecutive classes n 
and n + 1) that cannot be retrieved easily. In particular, they do not allow one to 
build the eigenvectors of a model knowing those of a simpler system. 

X We mention here that multispecies ASEPs with other local-interaction rules have been introduced, 
see [3[l0l[ll[28] 
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The above descriptions of the matrix ansatz and of the identification maps indicate 
that they operate in reverse directions and, hence, that they should be related. This 
also suggests that the matrix ansatz may have a range of applicability that exceeds 
stationary-state properties; the matrix ansatz should be useful also to build "excited 
states" of the A^-ASEP knowing those of the {N — 1)-ASEP. In this perspective, the 
matrix ansatz could be viewed as a 'lifting procedure' from the {N — 1)-ASEP to the 
A^-ASEP that creates a new species of particles by separating a given class into two 
consecutive subclasses. Further, we have seen that there are N different choices for 
the identification maps from the A/'-ASEP to the {N - 1)-ASEP. Are there also N 
different matrix product representations that would correspond to splitting a species 
n into two consecutive classes n and n + 1? Could there be N different generalizations 
of the Ferrari-Martin algorithm? 

The objective of the present work is to answer these questions in the TASEP case 
by reformulating them in the appropriate mathematical framework and stating them 
in rigorous terms. We will introduce a generalized matrix ansatz which constructs a 
lifting operator from the {N - 1)-TASEP to the A^-TASEP. This wih provide with a 
new and much broader perspective about this technique and shed light on the recursive 
structures that underlie the A'-TASEP dynamics. As a result, the problems stated 
above will be solved and families of quadratic algebras that encode these recursions 
will be constructed. The outline of this work is as follows. In section [2j we define 
the dynamical rules of the A^-ASEP, describe the characteristics and some spectral 
properties of the Markov matrix, define identification operators between different 
systems (sectors) and formulate the fundamental problem addressed here. From 
section [3] onwards, we restrict our consideration to the TASEP case, unless explicitly 
stated otherwise. In section |3| we recall the Ferrari-Martin algorithm of the stationary 
state of the A^-TASEP and describe the associated matrix product representation. 
The main results are derived in section [4] we construct a set of generalized matrix 
ansatz and we prove that this allows us to define a family of conjugation operators 
from the (A" — 1)-TASEP to the A'-TASEP; each of these matrix ansatz leads to a 
different quadratic algebra and we find explicit representations for all these algebras; 
finally, we find a generalized Ferrari-Martin algorithm corresponding to each algebra. 
Concluding remarks are given in section [5] The appendices contains examples and 
technical details. 



2. The A^-ASEP: Definition and properties 

2.1. Definition of the model 

The A'-ASEP on the ring Zl with L sites is defined by the following dynamical rules. 
Each site i G Z^, is assigned with a variable (local state) ki G {1, . . . , A" + 1} (A" > 0). 
We introduce a stochastic process such that nearest neighbor pairs of local states 
(a,/3) = {ki, ^i+i) are interchanged with the following transition rates: 

a3 ^ 3a < ' (2) 

^ ^ \ q if a> ^ ^ 

p and q being real nonnegative parameters, with the choice q < p without loss of 
generality [§J In particular, for g = the model is totally asymmetric and is called the 

§ In the paper 6 , the dynamical rule is defined as | ^ if a < ^' have to replace p ^ q or 

reverse the variables ki N -\- 2 — ki, if we wish to compare the present work with j6j . 
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7V-TASEP; for < g < p, the model is said partially asymmetric (7V-PASEP); q = p 
is the symmetric case (A/'-SSEP). 

The dynamics is formulated in terms of the continuous-time master equation for 
the probability P{ki • • - kL'^t) of finding a configuration ki ■ ■ ■ kL at time t: 



^P(/ci ■■■kL;t) = ^ 0{ki^i - ki)P{ki ■ • • ki-iki^ikiki^2 ■ ■ ■ t) 



dt 



Q{h-h+i)P(ki---kL]t), 



(3) 



where is a step function defined as 

r p (x<0), 

e(x) =1 (x = 0), (4) 
[ q (x>0). 

This model can be regarded as an interacting multispecies particle system by 
interpreting the local state A:^ = a as representing the site i occupied by a particle of 
the ath kind. In the conventional terminology, a particle of the ath kind for 1 < a < 
is called "ath-class particle." We can view the local state A' + 1 as a vacant site (or a 
hole), which is often denoted by 0. However, for later convenience, it is better to use 
A" + 1 instead of and we shall stick to that convention. 

2.2. The Markov matrix 

Let {|1),...,|A' + 1)} be the basis of the single-site space C^+^. The tensor product 
1^1 • • • ^l) = (8) • • • (8) ^ (^([^N+i^^^L corresponds to the configuration ki ■ ■ ■ kL 
on the ring. The probability vector at time t can be written as 

I^W)= E P{h---kL;t)\k,---kL). (5) 

l<ki<N-\-l 

In this language, the master equation (|3| becomes 

||P(t)) = MW|P(t)), (6) 

where the (total) Markov matrix M^^^ is of size (A" + 1)^ by (A" + 1)^. This Markov 
matrix can be written as the sum of local linear operators [M^^^j acting only 
on the ith and the {i + l)th components of the tensor product: 

V / 2,2 + 1 

ieZL 

where the action of the local Markov matrices on a bond is given by 

A^^oi = E - + ©(« - • (8) 

Q!,/3=l 

Note that off-diagonal elements of M^^^ are p, q or 0, and the diagonal elements 
are expressed as px + qy with nonpositive integers x and y. The sum of entries 
in each column of M^^^ is 0, assuring the conservation of the total probability 
^i<ki<N+i P{ki ■ - • kL;t). The master equation can be solved formally as 

|P(t))=e*^<"'|P(0)), (9) 
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and thus the eigenvalues and the right eigenvectors of M^^^ give information for 
physical properties of the model. The Markov matrix has the eigenvalue 0, and we 
call the corresponding right eigenvectors stationary states. All the other eigenvalues 
have strictly negative real parts (Perron- Frobenius theorem [45j), which characterize 
the relaxation to the stationary states. We remark that the Markov matrix M^^^ is 
Hermitian only for q = where it becomes the Hamiltonian of the 5/(A^ + l)-invariant 
Heisenberg spin chain. 



2.3. Particle conservation and Hasse diagram structure 

In view of the transition rule ([2|, the total Markov matrix obviously preserves the 
number of particles of each kind. Let > denote the number of particles of the 
ath class (respectively holes), ioi 1 < a < N (respectively a = N ^ 1): 

^•a2--^--- 7V + A/' + 1 . (10) 

m\ 7712 miv+1 

Then, the state space decomposes into sectors labeled by m = (mi, . . . ,mAr+i) with 
the constraint mi + • • • + mAr+i = L: 

(C^+i)«^ = 0^^. (11) 

m 

The Markov matrix has a block diagonal structure, leaving each sector invariant 

MW = Mm e EndKn, (12) 

m 

where the square matrix Mm acts on the vector space Vm spanned by all configurations 
belonging to the sector m: 

Kn= £\ki---kL). (13) 

#{i\ki=j]=mj 

The dimension of Vm is given by dim Vm = — r and the total number of sectors 

rn o J rn mi!---mAr+i' 

for given L and is ^ ~j\f^ 

Let us call a sector m = (mi, . . . , mAr+i) a basic sector if m^ > for all 
n = 1,...,A^ + 1 {i.e. there exists at least one particle of each type). Note that 
there are no basic sectors for N > L. Thus, we shall always take N < L — 1. For 
example, for L = 4, the list of all basic sectors for different values of N is given by 

(1,1,1,1), 
(2,1,1), (1,2,1), (1,1,2), 

(3,1), (2, 2), (1,3), ^'^^ 



N = 3 

N = 2 

N = 1 

N = 



(4). 



L-1 



The number of basic sectors for given L and N is given by ^ ^ J . Let A4 be the 

set of all labels for the basic sectors: 

M = {(mi, . . . ,mAr+i)|0 < TV < L - 1, m^ G N, mi H h mAr+i L}. (15) 

We now introduce an alternative labeling of the basic sectors [6 that will be very 
useful in the following: let sj be the total number of particles of classes k < i.e. 

Sj = mi + m2 + • • • + mj. (16) 
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We have rrij = Sj — Sj-i > with the convention = 0, and thus each basic sector 
of the A/'-ASEP can be labeled by the set 5 = {si, . . . ^sn} with < 5i < 52 < • • • < 
sn < L. The set 5 is a subset of 1^ = {1, 2, . . . , L — 1}, i.e. the set s is an element of 
5, the power set (the set of all subsets) of Q. We can identify A4 with S by the one- 
to-one correspondence (16). In the following, we shall use both labels equivalently: 
for instance, the invariant vector spaces (respectively the Markov matrices acting on 
them) will be denoted either by Vm or V5 (respectively or M5). 

An example of the identification M. ^ S for L = 4 is given in figure [1] The set S 
is equipped with a natural poset (partially ordered set) structure with respect to the 
inclusion C, which is encoded in the Hasse diagram [46 . In the present case, it is just 
the L — 1 dimensional hypercube, where each vertex corresponds to a sector. Every 
link of the hypercube becomes an arrow i ^ 5 meaning that t C 5 and #s = #t + 1. 
The "maximal sector" Q = {1,2,...,L — 1} corresponds to a unique sink and the 
"minimal sector" corresponds to a unique source, as in figure [TJ 




The following spectral properties in this Hasse diagram were proved in [6j: 

• Spectral inclusion: Let Spec (5) be the multiset of the eigenvalues of Mg with 
the multiplicity of an element representing the degree of its degeneracy. Then, 

Spec(5) D Spec(t) (17) 

for any pair of sectors 5 D t. In particular, Spec(r^) contains the eigenvalues of the 
Markov matrix Mg of all the sectors 5 G 5. An example of this spectral inclusion 
for L = 4 is displayed in figure [2j 

• Spectral duality: Using the spectral inclusion theorem, we can classify the 
eigenvalues E G Spec(5) into two types: (i) eigenvalues that already exist in 
lower sectors u C 5, (ii) eigenvalues that appear at 5. Let us call eigenvalues 
of the second type genuine eigenvalues. More explicitly, the genuine spectrum 
Spec° (5) of the sector 5 is defined as 

Spec°(0) := Spec(0) = {0}, Spec°(5) := Spec(5)\ |J Spec(u). (18) 

uCs 
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{1,2,3}: -4, (-3±i)2, (-3)^ -|±^, 
(-2)4, -|±A^,(-l±i)2, (-1)3,0 



{2,3}: -3±i , -3, (-2)^ 
_3_^ i V 
2 2 



^'iAil-.-lii, (-1)2,0 



> 






{3}: - 


2, -l±i,0 

> 



{1,3}: (-3)V|±>^ 



2 2 



(-2)2, (-l±i)2,-l,0 



{1,2}: -3±i , -3, (-2)", 

3 I i VT f 1 \2 n 

"2 2"' 




{2}: -3,-|±^ 

(-if, 




{1}: -2, -l±i,0 



0:0 



Figure 2. Spec (5) for L = 4 with p = 1 and g = (TASEP case). Each 
superscript denotes the multiphcity. For example, the Markov matrix Mj^^s} has 
eigenvalue —2, and its degree of degeneracy is 2. 



Then, the fohowing duahty relation was proved: 

Spec°(5) = + ^) - Spec°(5), (19) 

where the sector 5 = 1^ \ 5 is furthest from the sector 5 in the Hasse diagram. 

These two spectral relations were derived by introducing a family of linear 
operators connecting different sectors, that we now review. 

2.^. The identification operators cpis and the conjugation property 

Let us consider two (basic) sectors 5 = {si < - ■ - < sn} and t = 5\ {5^, . • . , Sn^}. We 
introduce a linear operator (pis : Vs ^ Vi, whose action on the basis vectors is given 
by 

• • • kL) eV^ ^ \k[--- k'^) e Vu with x' = x- #{i\ni < x}. (20) 

Note that (pis : Vs ^ Vi is surjective. For 5 = t, (pss is the identity operator for any 
sector 5. 



To understand the general definition (20), we give an example 5 = {2,3,5,8} D 



t = {2,5} with L = 9. These two sectors correspond to the following compositions of 
the system into different types of particles: 

2 3 5 8 

5:ll|2|33|444|5 

21 

t:ll|222|3333 ^ ^ 

2 5 
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According to these lists, we define (pis to be the operator replacing (or identifying) the 
local states according to the rules 3^2,4^3,5^3 (keeping 1 and 2 unchanged) 
within all the ket vectors • • • /c^) in For example, 

(/:)ts 1345214431) = |233213321). (22) 

More generally, for a pair of sectors t C 5, consider a chain 5o D 5i D • • • D 5^ of 
sectors such that 5o = 5,5^ = t and #5j = + 1 for all < j < u. Then, the 

decomposition 

^iS = V^tS^_i^5^_i5^_2 • • • ^S25i^Si5 (23) 

holds and is independent of the choice of the intermediate sectors [6]. 



2.5. The conjugation property 



A crucial property of the identification operators is that they provide a conjugation 
between different sectors. More precisely, for any two sectors such that 5 ^ t the 
following diagram commutes: 



Mi 



1 



Vi. (24) 



which means 



M^B = ^i.M, . (25) 
We shall call this type of relation a conjugation relation^ relating the dynamics 
governed by Mq and that governed by Mi. The identification operator (/^ts intertwines 
the two Markov matrices Mg and Mt, and we call ipi^ a conjugation matrix. The 
relation (25) was the key [6 in proving the spectral properties (17) and (19). Each 



right eigenvector \E) with eigenvalue E in the sector s can be projected down by LpiQ 
to the sector t, and (fi5\E) is also a right eigenvector with eigenvalue E in the sector 
5 (under the assumption ipi5\E) ^0). In particular, the stationary state of a sector 5 
is mapped by the identification operator to the stationary state of any sector t such 
that t C 5. 



2. 6. Looking for an inverse conjugation relation 



The identification operator (/^ts always maps an upper (more complex) sector 5 into a 
lower (simpler) sector t. This implies that identification matrix Lpi^ cannot help us to 
construct a right eigenvector in the sector 5 knowing the eigenvector with the same 
eigenvalue in a smaller sector t. It would be very useful if we could lift information 
from lower sectors to upper sectors in the Hasse diagram. One of the motivations of 
the present work can be formulated as follows: Can we construct a matrix ^/^^t which 
lifts up right eigenvectors from a lower sector t to an upper sector 5? This would 
be possible if we could define an operator ^jj^i from Vi to Vg such that the following 
diagram commutes: 



(26) 



Vi 



Mi 



> Vi 
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Equivalently, for any sectors t C 5 the conjugation matrix ip^t satisfies 

= M,^,i . (27) 
Note that the directions of the vertical arrows are opposite as compared to those in the 



previous diagram (24), which is a crucial difference. Equivalently, comparing equation 



(27) with equation (25), we observe that the order of the products is opposite. The 
property (27) has the following important consequence: Let \Ei) be any eigenvector 
of Mi with eigenvalue E. Then 

{^,i\E,)) = {M,\E,)) = E {^,i\E,)) , (28) 

which means that ?/^5t|^t) = is an eigenvector of Mq with the same eigenvalue E 
(if tJjsilEi) 7^ 0). In other words, the map ip^i allows to lift eigenvectors from a smaller 
sector to a larger sector and provides a constructive information from a simpler system 
to a more complex one (whereas the identification operator cpis erases information). 
The existence of tpsi is a very nontrivial property of the model. In the next section, 
we review the construction of the stationary state, which will be generalized later to 
define tpsi- 

3. The stationary state of the TV- TASEP 

From this section on, we set p = 1 and q = 0^ i.e. we consider the TASEP case, unless 
explicitly stated otherwise. In the A/"- TASEP dynamics, each particle can hop to its 
right nearest-neighbor site if the target site is empty or occupied by a higher-class 
particle. 

A stationary state is a right eigenvector of the Markov matrix corresponding to 
the eigenvalue 0. We denote it by \P) = P{r)\r): 

= mW|P). (29) 

For each sector m, the Markov matrix has a unique stationary state \Pm) up to 
an overall constant factor. For basic sectors 5 ^ m, we write IP5) = \Pm)- From 



the "grand-canonical" stationary state |P), which is a solution of equation (29), we 
extract the stationary state of each sector, by restricting the components of |P) to 
that sector. We emphasize that the multispecies TASEP on a ring does not satisfy the 
detailed-balance condition and exhibits non-vanishing currents in its stationary state, 
which is one of the major characteristics of nonequilibrium system^ 

The stationary state of the TV- TASEP is nontrivial as soon as > 2. 
The stationary state of the 2-TASEP was constructed by using a matrix product 
representation in [16]. However, this technique did not seem easily generalizable to 
higher values of N (see [37 for an attempt for the case N = 3). The solution to this 
problem came from two completely different directions. In [5 it was shown that the 
stationary state provided by the matrix product representation for the 2-TASEP could 
be interpreted in terms of weights in a binary tree. On the other hand, in [22 , these 
weights were rewritten in terms of Dyck paths which also appear as trajectories of a 
queueing process and therefore the 2-TASEP was reinterpreted as a queueing process. 
This fact was generalized to the A^-TASEP which was mapped into a system of coupled 
queueing processes [24 : this construction leads to the Ferrari-Martin algorithm for 
the stationary state of the A^- TASEP that we review in the next subsection. 

II However, it is important to keep in mind that boundary conditions are absolutely crucial in 
nonequilibrium physics: indeed, in a closed segment with reflecting boundaries, the detailed-balance 
condition is satisfied [TO] 
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3.1. Review of Ferrari and Martinis construction 

We now reformulate the algorithm found in [24 that constructs the stationary state 
of the A/"- TASEP, in terms of the convention adopted in the present work. The basic 
idea is to obtain the stationary state of 7V-TASEP from that of (TV - 1)-TASEP. 
The following algorithm is valid for any sector s = {si < ••• < sn}; In figure [sj 
we provide an explicit example for a ring of size L = 9 and with 5 = {2,5,6} {i.e. 
mi = 2, m2 = 3, ms = 1 and = 3). 




Figure 3. The Ferrari-Martin algorithm that constructs a configuration of A^- 
TASEP from that of (A^ - 1)-TASEP. This figure provides the specific example: 
F{wbbwbbbwb,231232313)=42U32U2. 

In order to relate a configuration of the {N — 1)-TASEP to a configuration of the 
A^-TASEP, we consider two lines, each of them corresponding to a lattice of size L. 

(i) On the upper line, we set sn black boxes b and (L — sn) white boxes w arbitrarily, 
that we write as ci • • • (with q = b^w). 

On the lower line, there is a configuration ki ■ ■ - ki, of the (A^ — 1)-TASEP 
corresponding to the sector s\ {sat}: thus, on the lower line, there are {s^ — s^-i) 
i^th-class particles (1 < u < N — I^sq = 0) and {L — sn-i) A'th-class particles 
{i.e. empty sites). 

(ii-1) Let {i^i \ . . . , ^si^} be the positions of the si first-class particles on the lower line. 
For the first first-class particle located at fi^id the nearest black box q/ = b 
with i^ < and put a particle of type 1 on it. If there is no such black box, 
put the particle 1 on the rightmost black box. For the second first-class particle 
located at the nearest unoccupied black box q/ = b with i' < 12'^ and 

put a particle of type 1 on it. If there is no such black box, put the particle of type 
1 on the rightmost unoccupied black box. Iterate this procedure Si times: i.e.. 
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find the nearest unoccupied black box = b witii i' < i\ for tiie ^tii first-class 

particle located at i^^\ and put a particle of type 1 on it or on the rightmost 
unoccupied black box if i' does not exist. 

(ii-2) Now we consider the second-class particles on the lower line. Recall that there 
are m2 = S2 — si of them, and set their positions as {zj^ , . . . , Zg^-sil- There are 
{sn — si) unoccupied black boxes remaining on the upper line. We must iterate 

the following procedure (52 — si) times: find the nearest unoccupied black box 

(2) 

Ci> = h with i' <i\ for the ^th second-class particle (1 < ^ < ^2 — Si), and put 
a particle of type 2 on it or on the rightmost unoccupied black box if i' does not 
exist. 

(ii-z/) In the same way, we go on for third-, fourth-, • • •, {N — l)th-class particles: for 
the z^th-class particles, there are {sj, — Sj,-i) z^th-class particles on the lower line, 
with positions \ . . . , 4l!-s^_i}- There are {sn — s^-i) unoccupied black boxes 
remaining on the upper line. Iterate the following procedure {sy — Sy-i) times: 
find the nearest unoccupied black box = h with i' < i[^^ for the ^th z^th-class 
particle {1 < i < Sj^ — s^y-i), and put a particle of type u on it or on the rightmost 
unoccupied black box if i^ does not exist. 

(iii) Now, there are {sn — sat-i) unoccupied black boxes that remain. Put particles 
of type N on them. 

(iv) Regarding the {L — sn) white boxes as an empty sites, i.e. as particles of type 
A/' + l, we have thus constructed a well-defined configuration F{ci • • • c^, ki ■ ■ ■ kL) 
of the A^- TASEP on the upper line, belonging to the sector 5, starting from a 
configuration ki • • • kL of the {N — 1)-TASEP on the lower line that was in the 
sector 5 \ {sat}. (Note that an N represents an empty site in the sector 5 \ {sat} 
whereas N -\- 1 does in 5.) 

Note that different configurations of {N — 1) -TASEP on the lower line can lead 
to a same configuration of the A/'-TASEP on the upper line: for example. 



It was proved in [24] that the stationary weight of a given configuration in the 
sector 5 is given (up to normalization) by the sum of the weights of all configurations 
in 5 \ {5 at} that are related to it through this construction. Equivalently, we have 



Here the summation ^ runs over ci ■ ■ ■ cl and ki ■ ■ ■ kL with = b} = sn and 

ki ■ ■ ■ kL belonging to the sector 5 \ {sat}. 

3.2. Matrix product representation for the stationary state 

The Ferrari-Martin algorithm was restated as a matrix product representation in [18], 
The basic idea of the matrix product representation is to express the stationary 
probability as the trace of a product of matrices over a suitable algebra. This 
technique, initially invented in [15 for the (one-species) ASEP with open boundaries, 
has been generalized to many stochastic interacting particle systems including discrete- 
time updates, a second-nearest neighbor interaction and non-conservative dynamics 



F{bbwwb, 32133) = F{bbwwb, 31323) = 21443. 



(30) 
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(see [To] for an exhaustive review). For the A^- TASEP on Z^,, the stationary weight 
of a configuration ji . . .Jl is given by the trace of an L-fold matrix product: 



PUl---jL)=TT 

or equivalently as 



X 



(AT) 



(AT)- 



(g)L 



|P) = Tr 



(32) 



(33) 



In [38], an exphcit solution for the operators Xj^^'s was found. It is given by the 
fohowing tensor product recursions 

N 



E 



iN,N) 
JK 



X 



(AT-l) 
K 



for 1 < J < A/- + 1 . 



(34) 



K=l 



with Y^Xl - ± ^^^2 
in the fohowing table: 



Y {i.e., X\^^ =X; 



^2 



1). The operators <^j^^^ are given 





1 ••• X-1 


X 


1 ••• X- 1 






X 





^0(iV-l) 


X + 1 




-(|0(Ar-i) 



(35) 



We read 5^ 



1 and 3^t 



for X < —2. The operators (5, e and A are 



the fundamental building blocks that are ubiquitous in the matrix ansatz technique 
[Toll . These three operators generate a quadratic algebra and satisfy the following 
relations: 

5e = 1, 5A = 0, Ae = 0. (36) 
A common representation of this algebra is given by the infinite dimensional matrices 



5 = 



( 



V 



1 

1 




7 



e = 



/ 

1 



V 



A = 



I 1 



V 







.(37) 



Remark: The Ferrari-Martin algorithm can not be easily defined for the PASEP 
because the directionality plays a crucial role in the algorithm. Nevertheless, the 
matrix product representation can readily be generalized to the PASEP case (p = 1 
and 7^ 0) as follows: in the table (35), we replace the operators (5, e and A by the ^q- 
deformed' operators (^g, and A^ that generate a quadratic algebra with the following 
relations 



^q^q - Q(^q^q = (l " Q)% ^q^q 



qAS^ 



q-> 



Aq^q 



qCqAq. 



(38) 



Then, the matrix product representation (33) with this deformation provides the 
stationary state of the X-PASEP, as was shown in [38j . 
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3. 3. Interpretation of the matrix ansatz as a linear mapping 



The recursion relation ( 34 ) implies that each stationary weight of the A/"- ASEP can be 
expressed as a linear combination of weights of the (A^ — 1)-ASEP [38]. More precisely, 
using equation we can write 

|PW) = *W^)|P(^-1)), (39) 

where the matrix ^(^'^) is defined in terms of its elements as 

Oi---J.|*(^-^)|fci---M=Tr[a(f,)...a5.f)^] (40) 

for 1 < ji < N -\- 1 and I < ki < N. Moreover the matrix ^(^'^) has the 
following property that we shall call sector specificity: suppose that the configuration 
ji-'-jL belongs to a basic sector 5 = {si < ••• < sn} and that the element 
Oi • • • JlI^*^^'^^!^! ' ' ' ^l) is nonzero, then the configuration ki ■ ■ ■ kL belongs to the 
sector 5 \ {sat}. Conversely, if ji ■ - ■ Jl belongs to 5 and ki - ■ ■ kL does not belong to 
5\{57v}, then (jl • • • JlI^^^'^^ 1^1 • • • ^l) = 0- We will prove this property in appendix 
B. This property allows us to consider the mapping i1^3^s\{sn}' ^s\{sn} ~^ defined 
as the restriction of ^(^'^) to the sectors s and s \ {sat}. This mapping provides us 
with a construction of the stationary state of the basic sector s by lifting up that of 
5\{sn}: 

\P.)=A,A{sn}\Pa{sn})' (41) 
Using this equation repeatedly, we have 

IPs) = '^S,5\{sn}^5\{sn},s\{sn-1,Sn} ' ' ' S 1 , ^2 } { S 1 } S 1 } I ^ ' ' ' (^2) 

where |1---1) is the only configuration of the 0-TASEP on a ring of size L (the 
configuration in which all sites are empty). Hence, the stationary state of the basic 
sector 5 = {5i,...,5Ar}is constructed along the way 

^ {Si} {Si, 52} ^ >5\ {sn-1,Sn} ^ 5 \ {sn} ^ 5 (43) 

in the Hasse diagram. 



On the other hand, using (23), we can project \Ps) down to the minimal sector 
via arbitrary intermediate sectors, i.e. for any {ni, . . . , un} = {1, . . . , N}, 

^05n-i'Pbn-ibn-2 ' ' ' ^b2b,^bib\Pb) = constant|l •••!), (44) 



where Sx = • • • , Sn^}. Comparing equations (42) and (44), we observe that the 



mappings play a role opposite to that of the (^'s. The ip^s are thus good candidates 



to be solutions to the conjugation relation (27). This property will be proved in the 
next section. 

We emphasize that the matrix ansatz that we have considered above allows us to 
construct the stationary state only along a very specific path in the Hasse diagram: 
in the sector 5 \ {^at}, there are mi = Si particles of class 1, m2 = ^2 — si particles 
of class 2,..., nfiN-i = ^at-i — sn-2 and nfiN = L — sn-i particles of class N. In the 
sector 5, the number rrij of particles of class j, with 1 < j < — 1 is the same, but 
there are new m'^ = sn — sn-i particles of class N and m^^-^i = L — sn particles of 
class A' + 1, with m'^j^^ + m'^ = ttin- Hence, when the sector 5 \ {sn} is lifted up to 
5, a new species is created by splitting the particles of class A" (that have the lowest 
priority) into two subspecies of class A" and A'+l. In the Ferrari-Martin algorithm, this 
means that a new species is created from the holes. However, in the Hasse diagram, 
there exist different paths between two connected but non- adjacent sectors. This 
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observation suggests that the matrix ansatz and the Ferrari-Martin algorithm should 
be generalized in order to construct the stationary state via arbitrary intermediate 
sectors. 

In the next section, we show that the matrix ansatz not only provides a way 
to write the stationary state but also allows one to define mappings that satisfy 
the conjugation relation ( [27| ) Therefore, this technique provides a tool to lift 
eigenvectors by intertwining the dynamics corresponding to different values of N. 
Besides, we shall also construct intertwining operators between arbitrary sectors, by 
defining and using more general quadratic algebras than those considered previously. 



4. Conjugation matrices from a generalized matrix ansatz 

In this section, we derive a general conjugation relation between A/'-TASEP models 
with different values of TV by using a generalized matrix ansatz which allows us 
to create a new class of of particles by splitting any intermediate species into two 
subspecies. 



4-1- Generalized quadratic algebra 

We define a family of rectangular matrices }i<^<7v ^^^^ + 1) x A^, 



indexed by the integer n. The elements <^j^^^ = {J\a''^'^^K) of the matrix a 
are operators given in the following table: 



JN,n) 



j\K 


1 ••• n-1 


n 


n + 1 ••• TV 


1 

n-1 


0e0l^(^-^-^) 


A^(^-i)0 


A^(^-i)0(50l^(^-^-') 
050l^(^-^) 


n 





^0(n-l)(g) 
-||(g)(Ar-n) 


050l^(^-^) 


n + 1 




1^(^-1)0 





n + 2 
7V + 1 




1^(^-3)0 


1^(^-2)0^01^(^-^-2) 



(45) 



where we read S<^ 



1 and (^01 



for X < —2. The fundamental 



operators 5, e and A satisfy the relations given in equation (36). In general, each 
element of a(^'"^) is either or an {N — l)-fold tensor product of 1,^4,5 or e. Some 
examples are given in appendix A. Note that for the case n = N ^ one retrieves the 
operators that were given in the table ([35]). 
Let us define 



(a(^'-))' 



(46) 



^ This interpretation of the matrix ansatz as an intertwining operator between different dynamics was 
already used in a recent study of a one-species TASEP with open boundaries and with annihilation 
[H]. There, conjugation matrices with respect to the system size were constructed, which allowed 
them to calculate the normalization factor and certain correlation functions (see 49 for a related 
approach) . 
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We shall show that the following relation is satisfied: 

In other words, for any value of n, the matrix ^(^'"') allows us to embed the system with 
{N — 1) classes of particles into the system with N classes of particles. Using ^(^'"'), 
we shall be able to construct sector specific conjugation operators that intertwine the 
dynamics between any two basic sectors along any path in the Hasse diagram. 

We now derive equation (47). The method used is an extension of the hat matrix 
technique, that was developed to prove various matrix product representations [26 , 10]. 
Suppose that, for each value of n, there exists an operator valued (A^ + 1) x matrix 
^(Ar,n) g^^]^ ^^^^ following identity, that we shall call the hat relation^ is satisfied 



(48) 



Then, equation ( |47| ) is a consequence of this relation. Indeed, from the relation (48), 

(Ar-i)> 



we obtain, taking the L-fold tensor product. 



leZL ieZL 
^(a(--))"'0a(--)^(a(--))^^'-'-^^ (49) 



ieZL 

= 0. 



ieZL 

Taking the trace of this relation on the space on which the operators a^^'"^^ act, 
and noting that local Markov matrices sum up to total Markov matrices, we obtain 
equation (47). We emphasize that the hat matrices a*^^'"^^ are used in the proof but 
do not appear in the final result (47). 



To summarize, the conjugation relation (47) follows from the hat relation (48) 
and in order to show that this latter relation exists, we need to specify the operators 
a^^'^). We claim that the hat relation (48) is satisfied for the choice 

a^^'^) = d^a^^'^) with dn = diag(l, • • • , 1, 0, • • • , 0). (50) 

n AT+l-n 

This explicit expression of a*^^'"^^ leads to closed quadratic relations between the 



elements ajK = {J\a^^^''^\K) (I < J < N ^ 1,1 < K < N) (Note that for simplicity, 
we write ajK instead of <^j^^^). Indeed, using the expressions of the local Markov 
matrices, we have 

^^^^ . . f -ajKaj'K' {J < J')^ 

{JJ'Wi^'J (aW^^^aW^)) \KK') = I {) ( J = J'), (51) 

{JJ'\ (aW-)0aW-)) Mi''^;'\KK') = | '^J^^^'^' + ajK^aj^K {K < ^^^^ 

Besides, using equation (50), we can calculate each element of the right hand side of 
(48| as 

f -CLJKCLJ'K' {J <n< J'), 

{Jf\ (a^^'^^^a^^'") - a^^'^^^a^^'")) \KK') = I ajKaj'K' {J > n > J^, (53) 

I (otherwise). 
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Substituting equations (51), (52) and (53) into the hat equation (48) leads to the 



following closed quadratic algebra generated by the operators ajx'- 





1{K < K') 


II {K > K') 


A{J<n<J') 






\ or n < J < J' J 






C ( J = J') 


CtJKCtJK' = CtjK'CtJK 




D {J>n>J') 


O^J'KCtJK' = CLJK'CLJ'K 




( J>J'>n \ 
\orn>J>J') 


C^JKCLJ'K' + CLj'KdJK' 
= ajK'CLj'K 


ClJ'KClJK' = 



(54) 



(Note that the cases A-II and C-II do not give any nontrivial relation, and E-II 
for K > K' is equivalent to B-I.) Showing that the matrix ^(^'"^)(46) satisfies the 



conjugation relation (|47|) therefore reduces to checking that <^jx^^'s defined by the 
table (45) actually give a representation for (54). Thus the proof of the conjugation 
relation (47) reduces to a purely mechanical procedure. We have checked this for 



several values of (A^, n) by using Mathematica. 



We emphasize that the key ingredient is the generalized hat relation (48) together 



with the ansatz ( 50 ) for the hat matrices a^^'^^ which allows one to define a quadratic 
algebra. The fact that the ansatz depends on the integer n, with 1 < n < A^, provides 
a family of quadratic algebras indexed by n. We also note that the conjugation relation 



(47) shows that the quadratic algebras defined in table (54) provide representations 
for the stationary state of the A^- TASEP: this is a much more concise (albeit more 
abstract) proof than the one given in [38]. 

4.2. Sector specificity 

We consider a basic sector 5 = {si, . . . , sat} and let C{5) be the set of all configurations 



of 5. We define %^^\(,^,y ^ K by 

for 1 < n,n' < A", ji ■ ■ ■ Jl ^ C{5) and ki ■ ■ ■ ^ C {5\ {sn'})- In fact, 
nothing but a sub-matrix of ^(^'^). Noting that M^Vs C K and M^v/^ ^xV^ 



(55) 



{N,n) 
s,s\{s^/} 



IS 



C 



deduce from equation (47) that 

which is the conjugation relation between 5 and 5 \ {s^/}. The following property 
implies that the conjugation matrix V^^^^j^ vanishes unless n' = n. More precisely, 
we have 

(i) If /ci • • • /cl (1 < ^£ < A^) does not belong to C (5 \ {s^}), then 



Tr 



0. (57) 



(ii) Equivalently, 



if Tr 



(Ar,n) 



^ 0, then ki---kL eC{5\ {sn}) . (58) 
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This statement will be proved in appendix B. Note that the sort sequences of the 
sectors 5 and 5 \ {s^} can be represented as 




•n + ln + 2--- N N^l---N^l, (59) 
n + 1 • • • TV - 1 ^ N • • • N ^ , (60) 

respectively. The index n specifies which kind of particles splits when the sector 
is lifted to 5. We write simply 



We recall that n appears explicitly in (50) for the hat matrix a^^'^\ leading to the 



quadratic algebra generated by the matrix elements of a'^^'^^ It is important to note 



that the sector specification property depends on the representation (45). 



We emphasize that the statements (57) or (58) do not guarantee that ips,s\{sri} 



is non- vanishing {i.e. meaningful). However, based on exact calculations for small 
system sizes using Mathematica, we shall conjecture the stronger property that the 
mapping is injective. 

4-3. Uniqueness. 

Conjugation operators T from V5\{5^} to Vs that satisfy M^T = TM^^i^^^^ are not 
unique. Indeed, because of the spectral inclusion Spec(5) D Spec(5 \ {s^}), there 
exist, in principle, at least diml^\|s^| such conjugation operators [8|. We conjecture, 
however, that the conjugation matrix '05,5\{s^} is unique if the following additional 
constraint is imposed. 

Uniqueness conjecture: The solution to Mg?/^^ = '0s,s\{s^}M5\|s^} is unique 

up to an overall constant factor when the following constraint is imposed: 

0'i---jL|V^s,s\W}|/^i---/^L) =0, (62) 

for all configurations ji • ■ ■ Jl ^ and ki - ■ ■ kL G C{$ \ {sn}) such that 3i such 

that ki -\- 1 < ji < n or n -\- 1 < ji < ki. 



The conjugation matrix ips,s\{sn} constructed by using the operators given in (45) 
satisfies this condition (see the example ( |A.3| )) 

4.4- Conjugation relation for s D t 

Now we turn to the construction of the conjugation matrix between Markov matrices 
of arbitrary sectors 5 and t such that 5 D i. Let us set 5 = {si < ••• < sn} and 



5 \ t = {s^^, . . . , Sn^}. Using the conjugation relation (56) between nearest-neighbor 
pairs repeatedly, we achieve the conjugation relation 

^,tMt = M,^,u (63) 

where 

with 5x = 5\ {sni , . . . , } = t U {sji^+i , • • • , } (for simplicity, we have omitted the 



superscripts in the ip mappings as in equation (61)). Each is constructed by 



Note that the uniqueness conjecture does not claim that the representation (|45| for the hat relation 
(|48|) is unique. 
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using with y = Ux — #{z\z < x^Uz < Ux}. A priori, this definition depends 

on the order chosen to enumerate the set 5 \ t. However, we conjecture that ip^i is 
independent from the choice of the path in the Hasse diagram. 

4.4- ^- Commutativity. The path independence in the Hasse diagram can be 
summarized by the fohowing statement, which we have checked for several values 
of (A^, n) by using Mathematica. 

Commutativity conjecture: For any sector $ = {si < • • • < sn} and for 

1 < X < y < N, 

i^5,5\{s^}i^5\{s^},5\{s^,Sy} = i^5,5\{Sy}i^B\{Sy},5\{s^,Sy}- (65) 

This conjecture implies that for two sectors t C 5, with 5 \ t = {s^, . . . , ^n^}, V^st 
is independent of the path in the Hasse diagram. In other words, for any reordering 
such that {ni, . . . , Uu} = {^1, • • • , tIu}^ 

with 5x=5\ {Sn^ , . . . , } and 5x =5\ {s^^ , . . . , 5^ J. 

Remark: In particular, we note that the stationary state IP^) of each sector 
5 = {si < ■ ■ - < Sn} can be obtained as 

\Pb) =^SSiV^3i32 •••V^SAr-23^,_iV^5^,_i0|l---l), (67) 

where we can chose the intermediate sectors Sx = arbitrarily. Since 

the stationary state is unique in each sector, the compositions with different sets of 
intermediate sectors must be the same up to a constant factor. This observation 
supports the commutativity conjecture. 

Example: We take 5 = {2, 3, 5, 7, 11, 13} and t = {2, 7, 13} so that 5 \ t = {3, 5, 11}. 
If one chooses 5^ = 11,^^2 = 5 and s^g = 3, then the intermediate conjugation 
matrices are constructed from a*^^'^\ a*^^'^^ and a*^^'^\ respectively, and this leads to 
the conjugation matrix i/j^t 

i^si = '0{2,3,5,7,11,13},{2,3,5,7,13}'0{2,3,5,7,13},{2,3,7,13}'0{2,3,7,13},{2,7,13}- (68) 

Taking 

— 5, — 3 and — 11, we obtain the following conjugation matrix: 

i^si = '0{2,3,5,7,11,13},{2,3,7,11,13}'0{2,3,7,11,13},{2,7,11,13}'0{2,7,11,13},{2,7,13}7 (69) 

where the intermediate conjugation matrices are constructed by a^^'^\ a*^^'^^ and a'^'^'^^ 
respectively. The commutativity conjecture implies that tpsi = 

4.5. Generalized Ferrari- Martin algorithm 

The algebraic construction that we have presented for all values of n, with 1 < n < A^, 
can be turned into an algorithm to calculate the stationary weights of the A^- TASEP, 
that generalizes the original Ferrari-Martin algorithm [24]. Figure [4] provides an 
example for {N^n) = (5,3). 

Given two sectors 5 = {si < ... < sn} and s \ {s^}, the matrix product 
representation of the conjugation matrix V^5,s\{s^} leads to the following generalization 
of the Ferrari-Marin algorithm. 
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(ii) 



^ — 7^ 1 — 7 ^ 



' 7^ — 



2U4 3) 5 1214 5 4 3 




(iii) 



2 V 5 y 5 V 1 k 6 1 y 5 6 y 2 




(iv) 



''6^^ y7yr)/?y?>(^'6';;/2 




Figure 4. The generalized Ferrari-Martin algorithm that constructs a 
configuration of A^-TASEP from that of (A^ - 1)-TASEP. This figure provides 
the example F(6k;i(;6w;^/;6k;6i(;6,24351214543)=25516415362. 



(i) Set Sn black boxes and {L — Sn) white boxes arbitrarily as ci • • • (q = 6, w) on 
the upper line, and a configuration ki - - - kL of the sector 5 \ {sn} on the lower 
line. There are 



{si, — Sjy-i) i^th-class particles (for 1 < u < n — 1^ 5o = 0), 
(5^+1 — Sn-i) nth-class particles, 

{siy-^i — Siy) i^th-class particles (for n -\- 1 < u < sn-\-i = 
on the lower line. 



(70) 
(71) 
(72) 



1 



As in the original Ferrari-Martin algorithm, put particles of z^th class from u 
to u = n — 1^ according to the rule (ii-^^). 

(ii-i/) There are {sn — s^-i) unoccupied black boxes. Let {ii, . . . ^is^-s^_i} be the 
positions of the z^th-class particles. Iterate the following procedure {s^, — Sj^-i) 
times; find the nearest unoccupied black box q/ = b with i' < ii for the ^th 
^(= ^ij, and put the particle v on it or on the rightmost unoccupied black box 
if does not exist. 

Put particles of {u + l)th class from u = Ntoiy = n-\-l, according to the rule 
{iii-u) which is opposite to (ii). 

{iii-u) There are (s^y+i — Sn) unoccupied white boxes. Let {zi, . . . , Zg^^-^-g^} be the 
positions of the z^th-class particles. Iterate the following procedure (s^^+i — 5^) 



times; find the nearest unoccupied white box Q/ 



with i' > ii for the ^th 



u{= ki^)^ and put the particle + 1 on it or on the leftmost unoccupied white box 
if i' does not exist. (Note that for u = we "put empty sites" (i.e., particles of 
(A/' + l)th class).) 
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(iv) There are {sn — Sn-i) unoccupied black boxes and (s^+i — Sn) unoccupied white 
boxes. Put n's and (n + l)'s on them, respectively. We have a configuration 
■ • • cl, ki ■ ■ ■ on the upper line. 

This generalized Ferrari- Mart in algorithm constructs the same conjugation 
matrix: 

= ^\F{ci"'CL,ki"'kL)){ki"'kLl (73) 

where ^ runs over ci • • • and ki • • - kL with = h} = Sn and ki ■ ■ ■ kL belonging 

to the sector 5 \ {sn}- 

This algorithm allows us to close the loop between this work and the previous 
articles [lH [181 [38] . these previous works, the original Ferrari-Martin algorithm, 
in which empty sites played a very special role, was used to construct a quadratic 
algebra to represent the stationary state. Here, we have found a family of matrix 
product representations that allow to split any given species into two subspecies (so 
that empty sites do not play a distinguished role anymore). 

To conclude this section, we precisely show the equivalence (73) between the 
matrix product representation and the generalized Ferrari-Martin algorithm (GFMA). 
We first note that each nonzero element of <^j^^^ (45) has the form 

= a^jV^ ■ ■ ■ • • • ®4S with a(^"J €{1,A, 5, e}. (74) 

(Note the shift in the subscripts that occurs for v > n). Let 

^=0C|m)>, {{m| = (0,..., 0,1,0,...), lM)> = {{Mr (75) 

/i>0 (M+l)th 

be the space on which the matrix ct^f^^J acts (thus afj^"'> acts on ^®(^-i)). The 
fundamental matrices A, 5 and e act on |//)) as 

The form (74) implies that its trace is again decomposed as 

f n ^(4:^:1— K.) (every 4-f)/0; 



Tr ('a^'^'"^---a^- 



(Ar,n)\ _ 



l<u<N 



(77) 



^at least one ~ ^ ) * 

Furthermore, we find 

1<1^<N 1<U<N n' >0 



l<u<N 



(79) 



(otherwise). 

The symbol * denotes the case where there exists (/ii, . . . /i^-i, Mn+i, • • • , Mat) such 
that 

^nCl ■ ■ ■ = afXl . ■ . aj.fr.! Jm'.)) = ^ ^.). (80) 

The set of numbers (/ii, . . . jin-i-, Mn+i, • • • , Mat) which satisfies this condition is unique 
if it exists. This uniqueness is true for basic sectors. On the other hand, we draw a 
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vertical line at each bond in the GFMA as in figure [5] The correspondence between the 
algorithm and the action of the operators can be understood by regarding as the 
number ja of arrows v ^ v {v < n) or u ^ u {v > n) crossing each vertical line: 5 
decreases the number /i of arrows, e increases /i, A tests whether /i = or not, and the 
identity operator 11 indeed does nothing. Then we find that, for given configurations 
ji ■ - ■ JL and ki - ■ -kL of A^-species and {N — l)-species sectors, the matrix product 
gives a unique pattern of arrows (if * is satisfied). This means that 



ji ' ' ' Jl can be obtained by using the GFMA from ki • • • ki,: 

F{ci---CL,ki---kL) =ji---jL with q = b{ji <n),w{ji >n + l). 



(81) 



This relation is true because the representation (45) obeys the rule of the GFMA. (For 
instance, the uth element of a[^u'^^ for ji.ki < n and < 7i — 1 is A. which means 
that ji can be put on the iih site of the upper line only when no arrow v ^ v crosses 
the vertical line between sites i and i + 1.) Thus, this graphical construction shows 
that the matrix product representation and the GFMA are euivalent. 

Figure [s] provides an example of the correspondence in the case (A^, n) = (3, 2): for 



jl • • • jio = 3211414433 and A:i • • • A:io = 1233321212, the actions of a^.^^^^^ and a 



^(3,2) 
ji ki ,3 



give trajectories of the numbers of arrows 1^1 and 3^4, respectively. Namely, the 
condition * is satisfied for these configurations: one can show 



(3,2) 



a 



• a 



• a 



(3,2) 
Jiofeio,! 
(3,2) 
Jiofeio,3 



|0)) = |0)), 



a 



(3,2) 

(3,2) 
'Jifei,3 ■ 



(3,2) 
■ *S'io/cio,l 
.(3,2) 



Ml) 



• • a 



'jiofcio,3 



= (m'i 7^ 1), 



Ims)) = (Ms 7^ 0). 



(82) 
(83) 



4 ^ ^/T)^' 4 V 4 




Figure 5. Correspondence between the representation ( |45| and the generahzed 
Ferrari-Martin algorithm. 



Recursive structures in the multispecies TASEP 



23 



5. Concluding remarks 

The multispecies exclusion process exhibits rich combinatorial properties that can be 
encoded in a recursive structure known as the Hasse diagram. The relevance of this 
diagram can be understood by using identification maps: if one blurs the difference 
between particles that belong to two adjacent classes, a simpler system is obtained 
that inherits the properties of the original system. Successive identifications allow one 
to reduce the initial multispecies system to the one-species system. In this procedure, 
information is irreversibly lost. 

In this paper, we have shown that the matrix product representation, which was 
originally used as a method for representing the stationary state, allows one to define 
conjugation operators that relate systems with different values of N in the TASEP 
case. Therefore, the matrix ansatz provides a method to lift information from a simple 
system to a more complex one: it allows one to calculate not only the stationary state 
but some excited states as well, that describe how the system relaxes towards its 
steady state. More precisely, each link in the Hasse diagram corresponds to a lifting 
operator and is associated to a different quadratic algebra. Hence, the A^- TASEP 
leads to families of algebras, connected through compatibility relations. We believe 
that this feature is general and that the investigation initiated here could be extended 
in the following directions: 

(i) The A/'-TASEP is an integrable model and its Markov matrix can be viewed 
one special member of a family of commuting transfer matrices for the Perk-Schultz 
model [39[|42]. One natural question is to study if the matrix product representation 
for such vertex models allows one to define conjugation operators. 

(ii) We also emphasize that we have studied here the A/'-TASEP only, where 
particles hop in one definite direction. It is natural to expect that our results could be 
extended to the A'-PASEP. A naive guess would be to start with the TASEP solution 



(45) and make the following replacements 

S^Sq, e^Cq, Aq (84) 



where Sq^Cq and Aq satisfy the g'-deformed quadratic relations (38). However, this 



guess is correct only for n = 1 and n = N but wrong for 1 < n < A'. The basic 



mathematical reason is that the g-deformed quadratic relations ( 38 ) are not stable by 
tensor product unless q = 0: di different approach seems to be required for solving the 
A'-PASEP in full generality. 
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Appendix A. An explicit example 



Here we write down the operator- valued rectangular matrix a'^^'^^ (45 ) for a few values 
of (AT, n). 





a^^'^^ = 2l A , a^^'^^ =20 A , (A.l) 



Recursive structures in the multispecies TASEP 

^(3,1) = ^(3,2) ^ ^(3,3) = 



24 



1 


^ l(g)l 


(5(g) 1 




1 


1(8)1 


^(g)l 


(5(g)^ ^ 


1 


/ 1(8)1 


(5(8)6 (5(8)1 ^ 


2 










2 





A0l 


A(g)(5 


2 





A<S)t A<S)S 


3 


e(g)A 







3 


e^A 


mA 





1 

3 








4 


\ 1(8)6 




mi J 


4 


\ e(g)e 


l(g)e 


l(g)l yi 


4 


\ 6(8)1 


1(8)6 1(8)1 / 



(A.2) 



^(6,3) = 



/ l-l-l-l-l 




e-l'A-A-A 
e-l-e-A-A 
e-l-l-e-A 

7 V 6-1-1-1-6 



^.6-l-l-l 
Al-l-l-l 


l-e-A-A-A 
1-e-e-A-A 
l-e-l-6-A 

1-6-1-1-6 



(5-1-1-1-1 

A-S-1-1-1 

A'A'1'1'1 

1-1-A-A-A 

I'l-e-A-A 

1-1-1-e-A 

1-1-1-1-6 



A-S-S-1-1 
A-A-5-1-1 


1'1'1'A'A 
1-1-S-e-A 

l-l-^-1.6 



A-S-M-1 
A-A-1-5-1 





1-1-1-1-A 
l.l.l.^-6 



(5-l-l-l-(5 \ 

A-S-l-M 

A'A'1'1'6 







l-l-l-l-l / 



(A.3) 



where we replaced (8) by We also give an example for the new conjugation matrix ip 
as well as the identification (p for {1,2,3} and {1,3} ((A/", n) = (3, 2)) with L = 4 (the 
entries of the matrices equal to are replaced by - for better readability): 

\ 



1 1 • 
1 • 

1 

1 



'0{1,2,3}{1,3} 



1 1 
• 1 



1 1 
1 1 



1 



' 1-1 



^{1,3}{1,2,3} 



1 1 



• 1 • 
111 
1 • 1 
1 • 
1 



1 

1 1 



1 1 
1 



(A.4) 



1/ 



(A.5) 
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where the bases are arranged as 1234,1243,. . .,4321 for {1, 2, 3}, and 1223,1232,. . .,3221 
for {1,3}. They satisfy 

V^{1,3}{1,2,3}^{1,2,3} = ^{1,3}^{1,3}{1,2,3}, (A.6) 
V^{1,2,3}{1,3}^{1,3} = ^{1,2,3}V^{1,2,3}{1,3}- (A.7) 

The Markov matrix M|i 3| has eigenvalue E = —1, and we write its corresponding 
eigenvector \E): 

\E) = |1223) - |1322) + |2132) - |2213) - |2231) + |2312) - |3122) + |3221) (A.8) 
The conjugation matrix '0{i,2,3}{i,3} lifts |^) to the sector {1,2,3} as 

1^0 '= V^{1,2,3}{1,3}|^) = 

|1234) - |1243) + 2|1324) + |1342) - 2|1423) - |1432) - |2134) + |2143) 

-2|2314) - |2341) + 2|2413) + |2431) + |3124) + 2|3142) - |3214) - 2|3241) ^^'^^ 

+ |3412) - |3421) - |4123) - 2|4132) + |4213) + 2|4231) - |4312) + |4321) 

The vector \E') is an eigenvector of M|i 2,3} (^{1,2,3} 1^0 ~ E\E^)^ . The identification 
operator reconstructs the eigenvector in the sector {1,3}: 

'P{l,3}{l,2,3}\E')=i\E)- (A.IO) 

Appendix B. Proof of the statement ( |58| ) 

We shall use the following property: let bi G {l,A,5,e} {1 < i < L). 

#{i\bi =S}j^ #{i\bi = e} ^ Tr (5i • • • 5l) = 0. (B.l) 



We also note the decomposition ( |77[ ). 

We first consider the case (A^, n) = (3, 2) as an example. Noting the 
properties (B.l) and (77), we find the following necessary condition such that 

Tr (ci^^^'^^ • • • 7^ 0- every af.'^^ is nonzero and 

for z/ = 1, 3. For given 3- and 2-species configurations ji ■ • ■ Jl and /ci • • • /c^, we have 

= 1} = # {^af^li = ^}+ = ki = 1}, (B.3) 

#m = 4} = # {i\af£s = e} + = h + l = 4}, (B.4) 

#{i\ki = 1} = # = e} + #m = h = 1}, (B.5) 

#{i\h = 3} = # {i\af^l:, = 5]+ - 1 = fci = 3}, (B.6) 

assuming that every a^^^' ^} 0. (See the explicit form for a*^^'^^ (A.2).) The condition 
(R2| and the relations' i|r3|-(|b]6| imply that 

= 1} = = 1}, W\h = 4} = j^{i% = 3}. (B.7) 

This consequence means that if ji • • • Jl ^ C{{si, 52, S3}) and Tr (ct^^^'^^ • • • 7^ 

then ki - ■ -kL G C({5i, S3}). 
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In the GFMA, the condition (B.2) corresponds to the conservation of the number 
of arrows when going around the ring (80). As in figure [sj for ji • • • Jl = 3211414433 
and ki ■ ■ ■ kL = 1233321212, we actually observe 

#oi4Si = ^} = #{^i4Si = ^} = 3> 



s} = 2. 



For the general case, a necessary condition for Tr (^<^j^'|^ 



(B.8) 
(B.9) 

^ is that 



every a^.^. is nonzero and 



= ^} = = 4 for V I. G {1, . . . , n - 1, n + 1, . . . , iV}. (B.IO) 



We can assume that every o!^^]^^^ is nonzero. From the definition (45), we find 
*{i\3i = v) = 




a^^.^ = e} + mji = h + l = v} {n + 2<v<N + l), 

= 4 + = ki = v} {l<v<n-l), 

= ^}+ -l = ki = p} in+l<p<N). 

Consequently, we have 

= v} = #{i\ki = v} for 1 < 1/ < n - 1, 
= v} = #{i\ki = 1/ - 1} forn + 2 < 1/ < + 1, 
which exactly means that if j\ - ■ ■ Jl ^ C (s) (s 




(B.ll) 



(B.12) 



Tr 



/ (JV,J7 



(B.13) 
(B.14) 

{51, . . . ,SAr}) and 
7^ 0, then ki ■ ■ ■ kL G s \ {sn}- (See the figures of the sort 



sequences (59) and (60)). We emphasize, however, that even if ji ■■■jL ^ C{5) and 



ki • • • kL ^ C {5\ {sn}), the matrix element Tr 
zero. 



{N,n) 



can still be equal to 
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